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We suggest that highly conducting oriented polymers with a fibril structure 
can be modeled by a regular lattice of disordered metallic wires with a ran- 
dom first-neighbor interwire coupling which mimics the cross-links between 
Q ■ fibrils. We show that such a model can be described by a non-linear sigma 

Q^ . model. Within a one-loop self-consistent approximation, we determine the 

position of the metal-insulator transition as a function of interwire cross-links 
concentration, interwire coupling and number of polymer chains in a wire. 

> 

\q ■ PACS numbers: 71.20Hk, 71.30.+h, 72.15.Rn 

m 

<N 
in 

s 

a 

o 
o 



13 



1 



Typeset using REVTeX 



Highly conducting doped polymers based on polyacetylene, polypyrrole, and polyaniline 
have recently attracted considerable interest. Due to the advent of improved chemical pro- 
cessing, their room-temperature can be comparable to that of copper. A T-independent 
Pauli susceptibility, a linear T-dependence of thermoelectric power and a large negative 
microwave dielectric constant further suggest that highly conducting polymers are intrinsi- 
cally metallic. However, in contrast to traditional metals, their conductivity decreases with 
temperature. This latter property is usually explained by the fact that highly conducting 
polymers are close to a metal-insulator transition (MIT) driven by disorder |l|]. The decrease 
of the conductivity with T then results from phonon-controlled localization effects. Moti- 
vated by these experimental results, many explanations of the transport properties of highly 
conducting polymers are based on dirty quasi-one- dimensional conductors (i.e. weakly cou- 
pled chain systems) [[|. The assumption of a periodic arrangement of the ID polymers 
chains which underlies such explanations is however quite unrealistic. Prigodin and Efetov 
(PE) have recently proposed a model which takes into account the irregular structure of 
the polymers ||. In highly conducting polymers, single chains are coupled into fibrils which 
are bent in space in a very complicated way. PE model the fibrils by weakly disordered 
metallic wires. The cross-links between the fibrils are described by interwire junctions. In 
the absence of junctions, all electronic states of the wires are localized by any weak disorder. 
PE have shown that the interwire junctions lead to a MIT and determined the position of 
the transition and the critical behavior. 

We suggest in this Letter that highly conducting oriented polymers with a fibril structure 
can be modeled by a regular lattice of disordered metallic wires with a random first-neighbor 
interwire coupling. The case where the wires contain a single chain has been recently studied 
numerically by Zambetaki, Economou and Evangelou (ZEE) [Q. We show that the MIT 
in this regular lattice of randomly coupled wires can be described by a non-linear sigma 
model (NLaM) which we solve in a one-loop self-consistent approximation. For a large 
number M of ID chains in the fibrils, we recover the main result of PE: an increase of the 
interwire coupling J favors a derealization of the electronic states. However we obtain a 
different critical behavior and the position of the transition differs from the result of PE. 
When M = 1, we show that a small value of J favors the metallic state while a large value 
of J induces localization in agreement with the numerical calculation of ZEE. 

We consider PE's model in the case of strongly oriented polymers and make the following 
Gedanken experiment: we stretch the wires in such a way that they form a regular (square) 
lattice of straight and parallel wires (Fig. [TJ). The interwire junctions now correspond to 
couplings between neighboring wires. In the case of a strongly oriented polymer, these 
couplings will have a short range. It is then possible, without inducing any qualitative 
change, to consider a simpler model where the couplings are allowed only between first- 
neighbor wires and in a direction which is perpendicular to the axes of the wires. To further 
specify the model, we assume that at each position x of a given wire, the coupling strength 
with a neighboring wire is equal to J with probability c (presence of an interwire junction) 
and vanishes with probability 1 — c (absence of a junction). For M — 1, this model exactly 
corresponds to the one studied numerically by ZEE. It is clear that the model is meaningful 
only if the probability to have a junction at a given position along a wire is weak (c 1). 
As pointed out in Ref . || , the existence of a delocalized phase in this model is a highly non- 
trivial phenomenon. Indeed, the random interwire coupling increases the dimensionality 
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(which favors the delocalized phase) but is also an additional source of scattering. 

We first consider the case where each fibril reduces to a single chain (M = 1). We will 
show below how the results can be easily extended to an arbitrary number of channels. We 
note r the elastic scattering time which results from intrachain disorder. In the absence of 
interchain coupling the electronic states in the chain are localized with a localization length 
Rq ~ I where I = vft is the mean free path and Vp the Fermi velocity along the chain. 
The mean-value of the coupling between neighboring chains is i± = cJ. Deviations from 
this mean- value are taken into account via a random hopping Vr\,(x) with second cumulant 
equal to 

Ki^V^x')) = t ± 5(x - x')5 hh±1 

where tj_ = ac(l — c)J 2 . a is the lattice spacing along the chains axis. I = (l y ,l z ) (l y ,l z 
integer) determines the position of a given chain. 8i u i 3 ±i means that l\ and I2 are first- 
neighbors. Higher order cumulants of the random variable V x are neglected. Averaging over 
the intrachain disorder and the random coupling V 1 - is performed by introducing N replica 
of the system. The partition function can be written as a functional integral over Grassmann 
variables ip a ^ (a = 1, ...,N). The action S contains two terms: i) Sq describing a weakly 
coupled chains system with a dispersion law ek = ufG^eI — kp) — 2i±(cos{k y d) + cos(k z d)) in 
presence of intrachain disorder (with elastic scattering time r) (hp is the Fermi momentum of 
a single chain, and d the spacing between chains); ii) S^ is coming from the random coupling 
V- 1 -. Using ([]]), we obtain 

SftB = — 7T Yl <V±i / dxdrdr' 

x 'W*(x, r)Vf (x, t')^,(x, t')W(x, r) 

-Vn*, r)4*(x, r')#(x, r)] , (2) 

where r,r' G [0, 1/T] are imaginary times. Eq. (0) shows that the averaging over V ± 
generates hoppings of particle-particle or particle-hole pairs between neighboring chains. 

In the following, we decouple the quartic term S^ is by means of a Hubbard- Stratonovich 
transformation and obtain the effective action of the low-energy fluctuations of the auxiliary 
field around its saddle point value. The parameters of this effective action are expressed as 
a function of t± and the correlation functions obtained from an action Sq which is related 
to Sq in a simple way (see Eqs. (|7|,|D below) 0. 

In order to take into account on the same footing the particle-particle and particle-hole 
channels, we introduce spinors 0, 4> defined by 
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and 4> = (C0) T where C is the charge conjugaison operator 0. ip^x) = ^\ 
Fourier transformed fields with respect to the imaginary time and u n = 7rT(2n +1). The 
action S^ m is rewritten as 
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SLfa <f>} = Jdx Tr[B(x)t ± B(x)] , (4) 

where we have introduced the matrix field B^ m (x) = 4>f n (x) <8> 4>i m (x). In (f|), tj_ should 
be understood as the matrix tu,v = ti3i,v±i (diagonal in the indices a,n,i where i = 1,2 
refers to the two components of the spinors). Tr denotes the trace over all discrete indices. 
Introducing an auxiliary matrix field Qfn m {x) to decouple the quartic term (|), we write the 
partition function as 

Z = j V<pV<p e - So[ ^ ] 

x JvQe~f dx b^iQ^^Q^+^W^Qi^]] ; (5) 

where tj 1 is the inverse matrix of t±. (One can verify that different decouplings of S^ is would 
not lead to any non-trivial phenomena.) The auxiliary field has the same structure as the 
field B(x) and therefore satisfies the conditions Q + = C T Q T C = Q 0. 

We first determine the value of Q in the saddle point approximation. Assuming a solution 
x) — Qo8 a ,/3fin,m3i,j we obtain the saddle point equations 

Qo = ^(q± = 0)-^EGf(k), 
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Gf (k) = iu n - e k + — sgn(cu n ) + 2iQ * , (6) 
L It j 

where £_i_(q_i_) = 2t±(cos(q y d) + cos(q z d)) is the Fourier transform of t±6i t ii±i. N]_ is the 
total number of chains and L the length of the chains. Eq. (|6|) is obtained using the Born 
approximation for the single-particle Green's function corresponding to Sq. From @, we 
obtain Q = (l/4r')sgn(a> n ) with 1/r' = 87riVi(0)t_i_ where iVi(O) = 1/ttvf is the density 
of states of a single chain. Within the saddle point approximation, we therefore obtain 
a change of the elastic scattering time due to the random interchain hopping. The total 
scattering rate becomes 1/r" = 1/r + 1/r'. 

We now consider the fluctuations around the saddle point solution. As in the standard 
localization problem ||, the low-energy fluctuations correspond to spatially slowly varying 
fields which satisfy the constraints Q 2 = Q sp2 and TrQ = 0. Following Ref. |J, we shift 
the field according to Q — > Q + Q SP — Q/2, with ijQf^ m (x) = S a ^S n ^ m 5ijuj n , and expand 
the action to lowest order in Q and Q. Using the saddle point conditions (^) we obtain the 
effective action of the Q field 



S[Q] = j dxTrlQ^t^Qix) - ttt ± l Q(x)] - J dx x dx 2 jiQi^ mn {x 1 ) ij n^(xil 1 ,x 2 l 2 ) ij Qff nm (x 2 ) 

(7) 

where a sum over repeated indices is implied. TZ^m * s defined by 



lZ a/3 (12)— ^1-^) ^) ^ftm(^i ^|1, -0 



\^(1,1|2,2) R^(1,2\2,1)J ' 

K,i(xih, X 2 1 2 \X 3 1 3 , X4Z4) = (^?U X ^hm(x2)^* Um (^ a i n (x3))s ■ (8) 
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The action So is equal to So + 2i / dx Tr[B(x)(Q SF — f2/2)]. For a system with time reversal 
symmetry, the Fourier transform 7£^(q) has all its components equal to i?"^(q). In the 
diffusive regime, it is easy to obtain the expression of i?"^(q) from the action So- In the 
replica limit (N — > 0), we have for the diffusive (Goldstone) modes {u) n u) m < 0) 

JCtoT'L = 2^(0) + + 8 ^ T "{ s ™ 2 (q y d/2) + sm 2 (g 2 rf/2))) . (9) 

The saddle point Qq introduces a "mass" term in the diffusive propagator. The exact 
propagator i?"^(q) can be expressed as in (H) but with renormalized diffusive coefficients 
D" and D'j_ replacing the bare coefficients v\t" and 8tj_T", respectively. Notice that the 
"mass" term 1/r' plays the same role as a finite frequency. R^(q) can be obtained either 
from So or from S ' but at a finite frequency 1/r'. Here S' Q ' is the action obtained from So 
by the replacement r — > r" . Eqs. (|7]f|) yield (after a rescaling of the fields in order to have 
Q 2 = 1 with 1 the unit matrix) 

S[Q] = l Nl (0)J2(D"q 2 x + (£>'[ + I) (sin 2 (^/2) + sin 2 M/2)))Tr[Q(q)Q(-q)] 

-|jV 1 (0)Tr[fiQ(q=0)]. (10) 

Notice that the "mass" term 1/r' in the propagator i?"^(q) is canceled by the uniform part 
of Tr[Q(x)t'J 1 Q(x)] (see Eq. (|7|)), i.e. the low-energy excitations correspond to diffusion 
(Goldstone) modes as expected from general symmetry arguments 0. The effective trans- 
verse diffusion coefficient appearing in fllQ|) is the sum of two terms. D'[_ oc t 2 is generated 
by single particle interchain hopping (fj_) and 1/r' oc i± comes from particle-particle or 
particle- hole pair-hopping as described by (g). The effective action fllCf), together with the 
constraints on the field Q, corresponds to a NLcrM for a system of weakly coupled chains. 
Expanding the sin 2 in (|T0|) to lowest order (which amounts to keeping only the lowest order 
gradient terms) would yield the standard effective action of an anisotropic disordered system 
0. Although this would be sufficient to obtain the critical behavior at the MIT, it is nec- 
essary to use the action (0) to determine the correct position of the transition. Anderson 
localization in weakly coupled chains systems has been studied by Prigodin and Firsov (PF) 
using a self-consistent diagrammatic treatment introduced by Vollhard and Wolfle 0. 
This latter approach is equivalent to a one-loop self consistent approximation of the effec- 
tive action fllpp . For a system of weakly coupled chains with interchain coupling t± e s and 
intrachain elastic scattering time r e ff, PF found that there is a MIT when 2t± e f[T e ff = 7 with 
7 a constant of order unity. The relation tj_ e ffr e g = const has been confirmed numerically 
T0fl . Since our final results will strongly depend on the precise value of 7 when M = 1, we 



will consider 7 as an adjustable parameter in the following. In the case of weakly coupled 
wires, the position of the MIT is determined by 2Mt± e f[T e s = 7 with M the number of chan- 
nels in the wires. This dependence on M can be obtained from standard scaling arguments 
nfl . Applying PF's results to the action (|T0|), we obtain that the position of the MIT is 
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The last step is to obtain the diffusion coefficients D" and D±. As mentioned above, 
these coefficients can be obtained from Sq but at a finite frequency 1/V. Again, we have to 
calculate diffusive coefficients in a weakly coupled chains (or wires) system so that we can 
use PF's results. The position of the MIT is determined by 2£_i_t" = 7 (2Mi±r" = 7 in the 
case of wires). In the localized phase {2i±r" < 7), the localization lengths are £ = -Rq/H an d 

= with I" = Vjpt". = \^2i±r"d and £j_ are the mean free path and localization 

length, respectively, in the transverse directions. In the case of wires, one should replace the 
one-dimensional localization length Rq ~ I" by the localization length of an isolated wire 
Rq ~ Ml". |e| = 1 — {2M/^)t±r" in the localized phase and vanishes in the metallic phase. 
The finite frequency diffusion coefficients are given by 

D"M = - VF \ ia ; DIM = - ^ • (12) 

1 + £ 2 Kk" 1 + e\^\r" 



D" and D'\_ are obtained from ( |I2"D with \u v \ = 1/t'. From (|II|l2D , we obtain the position 
of the MIT in the system of randomly coupled wires: 

r"/T> 8tjr" 2 2f_ n ^ 

1 + j£|V + / | e | V \ 2 M 2 ' { ' 

1 + MV (1 + ij2777 J 

Following Refs. H|J, we introduce the number of junctions per unit of length, p = 4c/ a, 
the number of junctions within the wire localization length, p = pRo, and the dimensionless 
interwire coupling a = (iraNi(0)J) 2 . For c C 1, ([Hp can be rewritten as 

1 A 7 2 , 2 



+ -7Z - = 



16a 2 M 2 

M 1 / M M \ 2 



• 4 = 1+ 2^ + M5l 1+ 2^-4WsJ ' (M) 

For M = 1, the position of the MIT strongly depends on the value of 7. For 7 < 1/2, 
p is a decreasing monotonous function of a. For 7 > 1/2, p decreases when a ^ 1 but 
increases when a ^ 1. Choosing 7 > 1/2, we reproduce the numerical results of ZEE 
which shows a non-monotonous behavior of p as a function of a (Fig. |2|). When a <C 1, 
r' > r on the transition line. The main source of scattering comes from the intrachain 
disorder. The main effect of an increase of a is then to increase the effective transverse 
diffusion coefficient D'[_ + 1/t', which favors the delocalized phase. When a. ^ 1, the random 
interchain coupling becomes the main source of scattering (r 7 < r). As shown in Fig. for 
7 > 1/2 this effect dominates over the increase of the effective transverse diffusion coefficient 
and favors the localized phase. We can also obtain from ( ^3[ ) the critical intrachain disorder 
W = (12t a ./r) 1 / 2 as a function of J (c fixed) or c (J fixed) (Fig. ^|). (The intrachain transfer 
integral t x is related to the Fermi velocity by Vp = 2t x a for a half-filled band.) The very 
good agreement between our results and the numerical results of ZEE (compare Figs. |^J^ 
with Ref. 0]) strongly supports the validity of our approach. 

For M> 1, the position of the MIT does not depend on the precise value of 7 (which 
should be of order unity). Contrary to the case M — 1, the main source of scattering 
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on the transition line remains the intrawire disorder even when a > 1. We find that p 
is a decreasing monotonous function of a (Fig. |2|). This result is in agreement with PE. 
However, the position of the transition does not correspond exactly to the one obtained by 
PE. In particular, we find that p approaches zero when a ^> 1 while PE obtained p ~ 1 in 
this limit. Moreover, since the transition is described by a NLcrM, we expect the critical 
behavior to be the usual one 0, in agreement with the numerical calculation of ZEE but in 
disagreement with PE's results. 

I am grateful to G. Montambaux and J. P. Pouget for useful discussions. I also wish 
to thank I. Zambetaki for sending me Ref. || prior to publication, C. Bourbonnais for 
discussions on the work reported in Ref. |5) and D. Boies for sending me a copy of his PhD 
thesis. 
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FIGURES 



FIG. 1. (a) Schematic representation of a random network of metallic disordered wires. Each 
intersection between solid lines represents an interwire junction, (b) Stretching the wires, we obtain 
a regular lattice. The interwire junctions become interwire couplings (dashed lines). 

FIG. 2. Phase diagram in the region c <C 1 for M = 1 (solid line), M = 10 and 20 (dashed 
lines) . 

FIG. 3. Critical intrachain disorder W = (12t x /r) 1 / 2 [t x = 1]. (a) vs J for c = 1, 0.5 and 
0.1 (from top to bottom), (b) vs c for J = 0.5, 3 and 10 (from left to right). Dashed-line: 
W = (24t x i ± /-f) 1 / 2 vs Jc. 
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